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I. INTRODUCTION 


Coherent states have proven to be a powerful tool in many areas of physics as well as mathematics. 
The name ‘coherent’ goes back to Glauber j^, who rediscovered Schrodinger’s states in the context of 
quantum optics. They are also used, for example, in geometric quantization j^, harmonic analysis and 
representation theory 0-0|. This broad applicability entailed a vast number of generalizations, just to 
mention some Bfl- 


In quantum gravity coherent states are employed to derive a semiclassical limit of the model in ques¬ 
tion. Especially in the absence of experimental data, this can provide important insights on quantization 


ambiguities and possible inconsistencies, 
so-called complexifier coherent states 


11 - 


n canonical loop quantum gravity (LQG) one uses, for example, 
151, going back to the pioneering work of Hall j^, in order to de¬ 


fine such a limit. For constraint systems such as gravity one has to decide on which space, the kinematical 
or the physical Hilbert space, the states shall be defined. Which strategy is chosen depends, of course, on 
the problem in question but in many cases it is easier to build coherent states on the kinematical rather 
than the physical Hilbert space. This truly applies to LQG where physical states are only known formally 
and brings in a new aspect that has to be respected as the states designed on the kinematical Hilbert 
space should not lose their properties when solving the constraints. 

From an heuristic point of view the implementation of a constraint H is related to a sort of time- 
evolution generated by H since 


V'phys‘=' =' J dt e'*^V’ 


( 1 ) 


gives a formal solution. In fact, many strategies such as group averaging and rigging map procedures (see 




f in order to solve 




models arose. 


e.g. |l6l |l take this as a starting point. Evidently, this ansatz is also advocating itse’ 
the Hamiltonian constraint of LQG and is the initial idea from which spin foam 
So instead of asking "Is a coherent state maintaining its coherence when solving the constraints?" it is 
tempting to simplify matters and ask: "Is the coherent state stable under the evolution generated by 
U{t) ■= or likewise "Is U(t)'ijjz still coherent?" These questions are as well of interest in quantum 

mechanics because mostly one is not only interested in the semiclassical behavior at a certain time but in 
the dynamical evolution. 

In this work, the necessary conditions for the existence of stable complexifier coherent states are inves¬ 
tigated. It is found that in general it is very hard to construct a complexifier adapted to the dynamics of 
a given model. Nevertheless, the derived criteria are form invariant under canonical transformation which 
opens the possibility to excess a wider class of models, namely, those that show a quasi-periodic motion. 


In the following section, the semiclassical properties of coherent states (section III HI) and the con- 

Thereafter, 


If 


struction principle of complexifier coherent states (section III Gl) are reviewed based on 
a stability criterium for finite dimensional models will be derived and in section IHI HI we will discuss a 
simplified ansatz to find solutions to this condition. Section [HI Gl contains a proof that it is, in fact, not 
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possible to use this simplified ansatz to determine other systems than the harmonic and radial oscillator 
that posses stable complexifier states. A generalized construction principle using so-called action-angle 
coordinates and the Hamilton-Jacobi approach is given in section flV Al and some examples are analyzed 
in section IIVBI The paper closes with a short discussion of the results. 

II. COHERENT STATES 

A. Preliminaries and conventions 

If not stated otherwise it will be assumed that the phase space AJ of a given system with finite number 
of degrees of freedom f is the cotangent bundle T*C of the configuration space C. The Hamiltonian vector 
field xj of a continuous differentiable function / on AJ is the vector field that satisfies the condition 
0 E where H is the symplectic 2-form on AJ and the Lie-derivative along Xf- The Poisson 

bracket corresponding to H is given by 

{/,£/} := ^{Xf,Xg) = Xfig] 

and multiple Poisson brackets are defined through the recursion relation {f, g}(n+i) •= {/){/) fl'}(n)} with 
{f,g}(^o) ■- g- The Liouville measure is the measure on AJ which is invariant under the action of the 
symplectic group that preserves H. 

Throughout the rest of this paper, x will denote the f-tuple (xi,•••,Xf), x-y = YJj^iXjUj the usual 
Euclidean scalar product and (p, q) a canonical conjugated pair, that is, they satisfy 

{Pj,qk} = Sjk and {pj,Pk} = 0 = {qj,qk} . 

Furthermore, O will denote a sub-algebra of the Poisson-algebra C°°(C) that separates the points of Al 
and T the complex conjugate. 

Under quantization we understand a map (Al,{',-},0) ->• {T-L, ^[■,-],0) where is a Hilbert space 
and O is a sub-algebra of the algebra of linear operators C{'H) on PL that is a representation of O. If 
not said otherwise the Hilbert space PL is the space of square integrable functions L^(C,d/i) on a suitable 
extension C of the configuration space with measure d/i. The scalar product on PL is usually given by 

= J dfi{x) f{x) g{x) . 


B. Semiclassical and coherent states 


In his lecture "Uber die Spektraltheorie der Elemente" |20|, given in 1920 at a meeting of the german 


physical society in Berlin, N. Bohr introduced the principle that the behavior of a quantum system should 
mimic the classical one for high energies. To formulate this statement in a more precise manner it is useful 
to introduce the notion of semiclassical states. This are elements hr the Hilbert space PL that are 
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associated to points m in A4. They are constructed in such a way that the expectation value of a quantum 
observable O in a given subalgebra O c C{'H) that separates the points of A4 is close to the classical value 
0(m) of the corresponding phase space function O. Stated differently, the following properties have to 
hold for all generic points m e A4 (i.e. points for which the denominator in ([2|),(|3|) and dH) is non zero): 


• Expectation value property 


• Ehrenfest property 


• Fluctuation property 


(V^tri) O tpm) 

0{m) 


ih{0,0'} 


{AuOlpm)'^ 


( 2 ) 


(3) 


(4) 


For most systems it is not possible to design semiclassical states for all observables simultaneously but 
it highly depends on the chosen subalgebra. A good example are the coherent states of the harmonic 
oscillator introduced by Schrodinger in 1926 2l| which have ‘good’ semiclassical properties for the linear 
span of the annihilator a, the creator and 1. In addition, these states have several other desirable 
properties which motivates the following definition. 


Definition 1 (Coherent states). A system of states is said to be coherent provided that in 

addition to and m the states also obey: 

• Overcompleteness (Resolution of identity): 

ln= dz/(m) |V’m)(V^m| (5) 

Jm 

for some measure v on AA. 


• Annihilation operator property: There exist operators z such that zifm = ■ 

• Minimal uncertainty: For the self-adjoint operators x = {z + z^)l2 and y = {z - z^)l{2i) the 
Heisenberg uncertainty relation is saturated, i.e. 

((i - = ((S - = ^|{[i,i])„.| (6) 


where (•),„ := (V’mrV’m)- 
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Peakedness property: For any Ai, the overlap function 




( 7 ) 


is concentrated in a phase space cell of Liouville volume 5l([j5, g])„ 


As stated in the introduction, coherent states have a broad application in many areas of physics and 
mathematics, which entailed a vast number of generalizations, so that by now ‘coherent state’ is not a 
clear-cut expression in the literature. In the subsequent section the generalization suggested in lll-ll5l| is 
reviewed, which in most cases preserves the properties mentioned in definition [1] 


C. The complexifier method 


The central point of semiclassical/coherent states is that they are continuously labeled by points in 
the classical phase space Ai. In the case of the harmonic oscillator, this is achieved by constructing the 
eigenstates of the operator a corresponding to the complex parametrization 


a := 


mu. i . 

— {q - p) 

I mu 


( 8 ) 


It is, of course, just one possible choice out of many different parametrizations of Ai and therefore many 
other state systems are imaginable that are coherent in the sense of definition [1] This is the starting point 
of the complexifier method introduced in lll-ll5l| that enables to directly relate a complex parametrization 
of Ai with a state system in T-L. 

Let C : A4 ^ M be a positive definite function with the dimension of an action that is smooth with 
respect to the Liouville measure on Ai, has a nowhere vanishing Hamiltonian vector field xc and grows 
stronger than linearly in p for each fixed point q e C. A function satisfying these properties is called a 
complexifier since 


Qi 


e g. 


(9) 


yields a complex coordinate system on Ai for a given parametrization {qi\i = of the configuration 

space C. The smoothness of C and the fact that xc is nowhere vanishing guarantee that C induces a 

non-degenerate, smooth transformation. If A4 = F*C then q^ defines a symplectomorphisrr0 that 

maps a point m e A4 to a point . 2 :(tn) in the complex extension C'^ of the configuration space. For example, 
2 

the complexifier C = generates the symplectomorphism 


Ai = 


C, iq,p) a = q-ip 


^ This is no longer true for M + T*C but Zi still provide good local coordinates due to Darboux’s theorem. 
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According to Bargman and Segal |J, l5[, the associated coherent states ijja also give rise to a transformation 
from L^(M, d/r) onto the space of square integrable, holomorphic functions on C through the 

integral transform 


f{q)^[Bf]{a) = / dfi{q)'tpa{q) f{q) 


( 10 ) 


If the measure v is fixed by the resolution of the identity ([5]) then (1101) is even unitary. Thus, the harmonic 

oscillator states play a similar role as plane waves for the Fourier transformation. In fact, for this specific 
2 

example with C = ^ one can show that up to a complex phase ijjaix) is equal to 




dfce"^ gifc(3;-a) 


( 11 ) 


Jj/^a 


where y ^ a denotes analytic continuation. 


Definition 2 (Complexifier coherent states). A coherent state associated to a complexifier C e C°°(A4) 
is an element of % of the form 


Mq) = [e 


where q' ^^(tn) denotes analytic continuation to 2 ;(m) = [e”*'^^c' g](m). 


( 12 ) 


This also explains why it was required that C has the dimension of an action and has to satisfy a growth 
condition. Namely, if Cfh would not be dimensionless then the exponential would not be well-defined. 
Apart from that, must decay fast enough to smooth out the divergence of 6. In equation (1111) this 

is achieved by the gaussian factor that is added through the action of . 

By construction, the states (|12l) are eigenstates of the annihilators 


Zi = e 




(13) 


resulting from quantizing (j9]) and therefore automatically obey minimal uncertainty ([6| . Also condition ([2]) 
is clearly satisfied for the subalgebra spanned by z and z. All the other properties required in definition [T]do 
not follow directly but are very plausible since the states are essentially a regularization of a distributions. 
Therefore, it is reasonable to expect that they generate a resolution of identity and are well-peaked on the 
phase space. 


III. STABILITY OF COHERENT STATES IN QUANTUM MECHANICS 

A. A stability criterion 

Recall that the classical evolution is generated by the flow of the Hamiltonian vector field XH of the 
Hamiltonian H. For a not explicitly time dependent Hamiltonian this means that a phase space function / 
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evolves as f{t) = e^* 


/(to) and the quantum evolution of such a model is described by the operator 


Definition 3 (Stable Coherent states). A coherent state '4’z{to) stable under the evolution t/(t,to) iff 


U{t,to)'4’z{to)=^"^^^^ Vti 


(14) 


where z{t) follows the classical motion of z{to) on the phase space. The system of coherent states, 
{V'z(m) 1^ ^called stable iff all states are stable. 


A simple example for a stable state system are the original coherent states fja of the harmonic oscillator. 
By a short calculation using ifa = ^ 1^) -^hoK) = (^ + 1/2)K) one finds 

Uho{t) fja = 'ifait) ■ 

Even more general, the set i s stable under the evolution generated by a HamiltoniaiJ^ that satisfies 


= if{d,t) (see e.g. 2J,l25|). For complexifier coherent states one can proof a very similar criterion: 


Theorem 1. Suppose the set of complexifier states ^to {'i/^(m(fo))l^(^o) ^ -M.} is over-complete and 
stable and the time evolution {/(tffo) is unitary then 

^%(to) = i/i(5i,---%,io) = (15) 

On the other hand, if equation (USD holds at to then there exist an €> 0 such that St is stable with respect 
to U{t, to) for all |to - t| < e. 


Proof. 

According to definition |3] the system Stg is stable iff 

Zj{to)U{t,to)\'lfz(to)) = Zj{t)U{t,to)\lfz(to)) 

^ [^%(io)]|V'z(to)) = ■ 

Here, Zj{t) := lJ\t,to)zj{to)lJ{t,to) is the time-dependent operator in the Heisenberg picture. Since •Sto 
is over-complete the last equation can only hold if ^Zj(to) is a function of the annihilators Zj{to), that is, 
if ^ Zj{to) = i fj{zi,... Zf,to) where i is just introduced for convenience. 

On the other hand, if equation (1151) holds then in some region |to ~ ^1 < ^ the Taylor expansion of z{t) 
exists and reads 

tn jn 

= E—377:%(^o) = Zj{to) + i {t - to) fj{zi{to),..., Zf{to),to) + ■■■ . 

T2/ • Clt' 

Therefore ij(t)|V’^(to)) = ^i(t)|V’^(to)) in that region. □ 


^ If H is required to be self-adjoint then / must be linear in a and H must be of the form H = ij{t)a^a + f{t)A + f* {t)a + j3{t). 
Note, there are even more general Hamiltonians under which certain proper subsets of coherent states are stable (see 

0,0). 






This implies the following two conditions that are necessary for coherent sates to be stable: 

Corollary 1. Stable coherent states can only exist if it is possible to find an complexifier C such that 
z = q obeys 

^Zj{t) = ifj{zi,...Zf,to) = (16) 

Corollary 2. If neither the complexifier C nor the Hamiltonian H are explicitly time-dependent then 
e^'^xc H must satisfy 

X ■ 

H = i ^jffjiqi, ...,qf)+ g{qi, ■■■,qf) (17) 

for some functions fj and g on C in order that the complexifier coherent state are stable. 

Proof. Since ^Zj{t) = {H,Zj} and since e~‘^^xc is a symplectomorphism it follows from corollary [1] that 

0 = {H,zj}-ifj{zi{t),...,zfit)) = H,qj} - i fj{q^... ,qf)] . 

Yet, {e*'^^c H,qj} is eqnal to fj{q^... ,qf) if and only if equation (1171) holds. □ 

To snmmarize: The existence of stable states is closely tied to the classical behavior of the model. The 
snbtlety hereby is that conditions m and (1171) are local. Even thongh it might be possible to find a 
complexifier satisfying these conditions in a neighborhood of a point {po,qo) e it does not ensure that 
the resnlting parametrization z = e~^^xc q is sensible. For example, z might not be defined everywhere on 
the phase space or it might be multivalued. Yet, for the system S to be over-complete it is important that 
z provides a ‘good’ parametrization: 

Lemma 1. If the system S is over complete then {zj} is a proper parametrization of Ai in the sense that 
almost every point m e A4 is uniquely determined by a fituple (zi,---,Zf). Here, ‘almost’ means up to a set 
of measure zero with respect to the measure ©. 

Proof. If this were not the case then there would exist a subset U in Ai with a parametrization w such 
that states V’w(m) associated to m e t/ cannot be expanded in terms of states in S. But then 

/ dl/(m) |V'z(m))(l/’z(m)l ^ Iw 

Jm-U 

which contradicts the assumptions. □ 

Note, this does not exclude parametrizations that break down at single points nor does it exclude 
parametrization that differ in region A from this in region B if the resolution can be modified accordingly, 
that is, if splits into + f^. Especially the last point is of interest since in many models the dynamics 
can be fundamentally different in separated regions and therefore it might be necessary to consider different 
complexifiers for those regions. It is even possible that in certain areas of Ai one cannot find stable states 
while for others there exist some. For example, inside a finite potential well the motion is periodic which 
suggest that there exist stable states while outside the well the motion is unbounded. This issues will play 
a significant role in section IIVI For now the main aim is to find solutions to (|17l) . 
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B. First solutions 


For a one-dimensional not-explicitly time dependent system equation (1171) reduces to 


E = ipf{Q)+9{q) , 

«_n 


(18) 


72=0 


where {C,= {C, {C, and / and g are arbitrary smooth functions. To find a first solution 

of this equation consider the following ansatz: 


H=^ + V{q), C=^ + U{q) 

2m 2muj 

and {C^,-ff}( 2 iv+i) = 0 for some > 1 


(19) 

( 20 ) 


Here, m is the mass, V and U are twice continuously differentiable real valued functions and w > 0 is a 
free parameter. Note that the even (odd) Poisson brackets in (I19p can only contain summands of even 
(odd) powers of the momentum. Therefore, / and g in (1191) have to be real and satisfy 

E7^{C.if)c2„)=9(9). (22) 

n=0 

A model with f degrees of freedom can allow at most f independent commuting Hamiltonian vector fields 
(first integrals of motion) which is why for f = 1 condition (1201) is equivalent to 

N 

{C,27}(27v) = E/3iCl^' • (23) 

j=0 

The restriction on the degree of the polynomial is due to the fact that {C, H}^^ 2 N) can be at most a 
polynomial of degree 2N in p. For = 1 equations (1211) , (|22l) and (1231) are replaced by 


{C,H}=pf{q), 

(24) 


(25) 

and {C, F7}(2) = l3iC + f3o . 

(26) 


For this choice the first two conditions in (1191) are equivalent to 

{uj-^V'-U')/m = f and ^ f - f U' = 2{H - g) 

2muj 

where ' denotes the derivative with respect to q. Thus f must be equal to a; and /3i must equal 2u. 
Inserting this into (|26p yields 


f U'=-2u:[U + h] 
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which is solved by f = toq + a and U = ^{q + a/co) ^ + (^2 for some real constants A and a. Since a just 
dehnes a shift in the conhguration variable it can be set to zero without loos of generality. Concluding, 

.,2 


= 


V 

2m 


TTl 2 2 \ -2 

+ —ojq +Xq 


C 


^ -2 

-+ —q 

2mL0 oj 


P 


(27) 


solves dm and (1201) for iV = 1. If A is zero this reduces to the usual Hamiltonian and complexifier of the 
Harmonic oscillator. The term X/q^ can be interpreted as the angular momentum contribution of a two 
dimensional oscillator with constant radius which justihes the name radial oscillator. The complexihcation 


z := e q can be computed by solving the ‘equation of motion’ 


ds 


q={C,q} 


for g(s) := e^'^^c q and then extending it analytically (s ^ -i). Instead of integrating equation 
can integrate 


(28) 


one 


ds 


2 Ids/ 


mu 


^ -2 


2 Us/ 


mu 


^ -2 

— 9 


c 

mu 


since gj C* = 0 and {C,q} = pl{mu). This hnally gives 


z = qis = -i) 




-2 

— 9 


-.1/2 


mu 


Despite that z depends on a square root, one can derive a well-dehned quantum operator whose eigen¬ 
functions 5y\y^z can be expressed in terms of modihed Bessel functions of the second kind. 

Furthermore, the obtained states share all desired properties, that is, they form an over-complete set, 
minimize the uncertainty, have small fluctuations and are of course stable under the dynamics. More 
details on the quantization of this system can be found in appendix 

C. A no-go theorem 


After the method applied in the preceding section has been proven so successful it is worth testing 
whether more solutions to (|T8l) can be found by generalizing the ansatz (I19p . Given two real valued, 
strictly positive, continuously differentiable functions a and /? on the conhguration space let 


H = a{q)^ + V{q) and C =]^[l3{q)pf+ U{q) 


(29) 


and suppose {C*,/7}(27v+i) = 0 for some N > 1 . 

Since (3 should be non-zero in order that the complexiher is well-dehned there exist a canonical transfor¬ 
mation sending 

rq 

p-^pf3{q) and dx [/3(3:)]” 

Jo 
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so that the first line in (1291) may be exchanged through 

„2 

H = a{q) + V (q) and C= — + U{q) (30) 

without loos of generality. To compute the multiple Poisson brackets of fl30|) it turns out to be handy to 
introduce the operators X = and P = -U'^ which yields 

n-1 n-2 

HPX'" H + Y. PX'^ PX^^ H + ■■■ . (31) 

u=0 fi+u=0 

The general strategy in the subsequent analysis is to reorder the terms by powers in the momentum. The 
degree in p of the terms PX'^^--- only depends on the number of operators jjO that are applied, i.e. 

degpX'^^Px''^-V = ^X-^P 

and 

degp X’'^ Px’^^-ap'^ = ^X-^P + 2 

where degp denotes the degree in p. Note that the total number of operators, |j X + [j P, corresponds to 
the grade n of the bracket {C, For this reason the bracket either contains only even or odd powers 

in p depending on whether n is even or odd respectively. This can be made more explicit by replacing 
(f3T]) through 




/i€N 

0<(n+2)/2 


The coefficients fn+ 2 - 2 n derived from (|3ip . that is. 


n-1 / o s.n-l-u 


(32) 


(33) 


where is the nth derivative of a function g with respect to q and the upper index of revers to 
the total number of operators while the lower one indicates the power in p. Instead of trying to directly 
compute these coefficients it is much more useful to consider the following recursion relation: 

/m = ^/m-i-t^'(2m + l)/”;\, m<n + 2, (34) 

with /q E 1/ and = a. 

Furthermore, due to | degp {C,H}(^ 2 N) < N + 1 and {C, H}(^ 2 N+i) = 0 finds that 


Af+l 

{C,H}^2N)= E amC^ 

m=0 


( 35 ) 


12 


for some constants am- Apart from that, all summands in (1181) whose degree in p is exceeding one must 
cancel which implies that 


N 


E 

n=m-l 


^/ 2 t = 0 and 
(in)! 


7V-1 f 

y — 

n=m-l (2n + l)! 



(36) 


for all m > 1. Note, there can be only one non-trivial term of power 2N + 2 and one of 2N + 1 which is 
why /|jv +2 /liv+i have to vanish. This in turn implies deg^ {C, ^ ^ deg^ a < 2N -2 since 

/liV+i proportional to The next non-trivial contributions are those of degree 2N and 2N - 1. 

Due to equation (1361) the coefficients and /Ijv-i ™ust satisfy 


/|^ = (2iV)(2iV-l)/|^-2 and fl^ll = i2N - l)i2N - 2)fl^:f . (37) 


Yet, taking the derivative of the equation on the right hand side and inserting equality (|34p yields 

^ = (2JV - 1)(2N - 2)^ /l^:^ ^ /l^ = (2N - 1)(2JV - 2)/|^-2 

which obviously contradicts the first equation in (|37p . Thus, must vanish for all m. With the next 
term one can proceed in the same manner. Combining (1361) . (1341) and = 0 gives 

n N f \n N-1 { \n 

n _ ^ ^ V ) f2n _ V' V / f2n+l 

“ dq n^-2 (2n)!^27V-2 - ^ (2n)!^27V-i • 

This in turn contradicts the second equation in (|37l) and therefore implies f^_i = 0 for all m which also 
reduces the degree of the polynomials a and {H, C} 2 n and proves 


Lemma 2. For N > 1 there exist an integer M with 2 < M < 2N - 1 such that fm = 0 for all m> M and 
all n. 


The above reasoning can be repeated to show that also summands of lower degree in the momentum have 
to vanish. But note, the number of conditions needed to derive a contradiction is increasing when the degree 
in p is decreasing since there are more and more non-trivial coefficients that contribute. Nevertheless, one 
can obtain a new condition for the coefficients from 


(-Y 

n=m-2 \^^)- 

by taking the derivative and using (I34p . More specifically, this yields 

N-1 ^ — 


0= E 

n=m- 

N-] 

= E 


n=m -2 (2rr)! 
N-1 ( — 


(2n)! 


[/lr/i-(2m + 2)[/72tr2] 


/ir/i Vm>0. 


n=m-2 

From this one can deduce another condition 


(38) 


(39) 


N-1 


E 

n=m-2 


i-r 

(2n)! 


n2n+2 

J2m+2 


V/i>0 


(40) 
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by first taking the derivative and then applying (1341) and (1391) and so forth. Repeating this procednre 
also for the odd terms generates 2m - 1 independent conditions for the terms of power 2m and 2m + 1. 
On the other hand, the nnmber of non-trivial terms |j in a tower of constant power m increases as 
jt - m + 1 and jj f^m+i = N-m + 2. Since the procednre breaks np as soon as [j is greater than 

the nnmber of available conditions, namely if It/ 2 m/ 2 m+i - 2m- 1, another trick is needed to eliminate 
more coefficients. 

2 rn 

Lemma 3. Let (7 = ^ + U{q), F = Y, FnC'^ and let G = G{p,q) be a phase space function that is 

n=0 

polynomial in p. 

(a.) Ifdegp G = 2m + 1 then {G,G} = F{G) has a solution iffU = Ciq + CQ for some constants ci,co e M. 

(b.) If degp G = 2m - 1 then {G^G} = F{G) has a solution iff U = ci{q + 02 )”^ + cq for some constants 
C2, Cl, Co e M. 


Proof. 

Proof of (a.): Recall that the complexifier is the only constant of motion with respect to the flow generated 
by G itself so that the homogeneous solutions to the first order PDE {G,G} = F{G) are functions of G. 
Therefore, the most general solution G with degp G - 2m + 1 is of the form 

G = F{G)g{p,q)+gh{G) . 


Here, gh is any polynomial of G whose degree is not exceeding deg^ F and g = g(p, q) is a linear function 
in p for which 


1 dg dg , 
{C,g} =Pj--^U . 
oq op 


Since U and are independent of p this can be only solved if g and U' are constant in q. 

Proof of (b.): Without loos of generality one can assume that Fq = 0 because the Poisson bracket does not 
change when C is shifted by a constant. Suppose g = phi(q) + ho^q) then 


G = 


nc) 

G 


9{p,q)+9h{C) 


is a generic solution that has the right degree in p iff 


C={G,g}=p{ph[ + h'o)-hiU' . 


This implies that /iq is constant, hi = \{x + C 2 ) and 2U = -{x + C 2 )U' which is solved hy U = ci {q + C 2 ) ^ + 
Co- □ 


Remember that lemma[2]guarantees the existence of a non-zero integer mo < N, for which deg^ {F, II}( 2 N) 
mo and degp {G, II}(^ 2 N-i) = 2mo ± 1, and that the potential U must be either linear in q (for 2mo + 1) 
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or proportional to g ^ (for 2mo - 1) according to lemma |3l The degree in q of the functions //f can be 
directly determined from equation (|33p . that is, 

degg fn +2 = degg a-n 


and 


degg fn+ 2 - 2 m = “ax [degg a-n + m + m deg^ U', deg^ ^ - re + m - 1 + (m - 1) deg^ U'] . 

Suppose the coefficients /2mo+i are not zero then they have to be constant as f 2 mo +2 ~ d. In this case, 
lemma |3] states that also U' has to be constant which is why 

0 = degg 72 ^ 0+1 = degg a - (2rero - 1) 


and 


0 = degq = max [deg^ a - (2rero + 1) + 1, deg^ V - (2mo + 1) + 1 - l] . 

This forces deg^ a = 2mo - 1 and deg^ V = 2mo + 1. Yet, 

0 = degg = max [2rreo - 1 - (2rreo + 3) + 2, 2mo - (2rreo + 3) + 1] < 0 

leads to an inconsistency and thus /Imo+i has to vanish which means that must be constant. By 
using lemma |3](b.) one can deduce that the potential U is proportional to q~‘^ and consequently deg^ a 
must be equal to 2rreo - 2 and deg^ V must equal 2mQ. Again 

0 = < ieg ,/“»*"<0 

leads to a contradiction so that all f^mo have to vanish. By repeating this argument one can finally show 
that must vanish for all rre > 2. Hence, is constant and = -f^x + for some e M 
which can be only achieved if a is constant as well and U is proportional to x~‘^. Obviously, the brackets 
{ff, C}( 2 n) s-rid {H,C}^^ 2 m)i can only differ by an over-all constant so that already {iif, C}( 3 ) = 0. 

This proves the following theorem. 

Theorem 2. Suppose the Hamiltonian and the complexifier are quadratic in the momentum p then the only 
system that solves {C, H}(^ 2 N+i) = 0 for some > 0 and (jlSp is the radial oscillator (I27D and canonical 
conjugates thereof. Furthermore, N equals 1. 

Note, all the equations used to prove this theorem are directly related to Poisson brackets and for this 
reason only hold up to canonical transformations. This point will be exploited heavily in the next section 
to derive a more general construction principle for stable coherent state system. 
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IV. ADAPTED COMPLEXIFIERS FOR INTEGRABLE SYSTEMS 

The preceding investigations have revealed that it is in general very hard to construct a complexifier 
adapted to the dynamics of a given model. Nevertheless, the derived criteria are form invariant under 
canonical transformation which opens the possibility to excess a wider class of models than those examined 
above. A common feature of the harmonic and the radial oscillator is their periodic motion which is why 
integrable systems that show a quasi-periodic motion seem to be especially promising candidates. This idea 
is to be elucidated in more detail in the subsequent section. In the first part a generalized construction 
principle for adapted complexifiers will be worked out using so-called action-angle coordinates and the 
Hamilton-Jacobi approach. This will then be tested on several examples. 


A. Generalized construction principle 


Throughout this section it will be assumed that the mechanical models in question are not explicitly 
time-depended and posses only a finite number f of degrees of freedom. 

A first integral of motioi^ is a function / on the phase space Ai that Poisson commutes with the 
Hamiltonian H, i.e {H,f} = 0 on the entire phase space. Two functions g,f € C^{Ai) are said to be in 
involution if {g,f} = 0 and (functionally) independent on a subset [/ c AJ if the one-forms df and d^ are 
linearly independent on U. If / and g are independent then in particular df a dg is non-zero on U. 


Definition 4 (Integrable System). A system with f degrees of freedom is integrable if there exist f first 
integrals Hj, j = in involution that are independent on a dense subset ofAi. 


The name is motivated by the fact that such systems are integrable by quadratures, that is, they are 
solvable by a finite number of algebraic operations. This insight goes back to Liouville and was later 
enlarged by Arnold by the following (see e.g. 2^): If the system is integrable then the level sets 


Alh = {nr e Ai\Hj{m) = hj} 


(41) 


are smooth submanifold that are invariant under the phase flow generated by the Hamiltonian H. If Mh 
is compact and connected then it is diffeomorphic to the torus pf and the motion is conditionally periodic. 
This means that Mh has coordinates Qj which parametrize the circles in and which evolve as 

^=u.(h). (42) 

For this statement to hold it actually suffices that the first integrals are independent on Mh. Furthermore, 
this angle coordinates can be used to parametrize the phase space (in the neighborhood of the invariant 


^ The term ‘first integral’ often refers to a global property while ‘constant of motion’ is used in a more local context. Yet, 
the nomenclature is far from being unique; here the conventions of will be used. 
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torus). Its conjugated momenta Ij can be found by a canonical transformation that, by definition, leaves 
the symplectic structure invariant. So 

f f 

^ dpj A dqj = ^ dij A d0j . 
i=i 

The parameters (I, 0) are called action-angle coordinates in the literature and are widely used in classical 
perturbation theory (see e.g. 2a, [27| for more details). Here, their simple time-dependence is of interest 
which is given by 


^ = 0 and — = cu(I) . 
dt dt 

This immediately shows that the complex parametrization, 

Wj := ^yij e*^-’ and Wj ;= \/Ij 
is ‘stable’ in the sense that it obeys 


dw,- 


= {H,Wj} = iiOj{l)wj . 


In addition, the pair (w, w) is canonical conjugated as 

{wj,Wk} = iSjk and {wj, w^} = {wj, w^} = 0 
and the polar decomposition of Wj reads 


(43) 


(44) 


(45) 




where 


Qj = \/2Ij cosOj and Pj = -\j2Ij sinOj . 

Since also (P,Q) are canonical conjugated Wj can be written as a complexifier coordinate. 


1 


w. 


j - Qj , 

V2 


with C = • P. But note, the parametrization given by (|44p does not satisfy the criteria of theorem 

[1] since the frequencie^ ojj still depend on Ij = WjWj. If the system is non-degenerate, which is the 
case if det I 0, then the invariant tori are uniquely defined and independent of the initial choice of 
the coordinates ( 0 , 1 ). This implies that, no matter which action-angle coordinates are selected, the 
frequencies will always depend on I. In other words (1441) can only give rise to stable states if the system is 
degenerate. More specifically, if H = for constant ojj then the associated coherent states are stable 


du 


because = 0 for all j, k. 


Tn trptipral- ^ 

in general, g^^ g^_ g^^ . 
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Before bothering about degeneracy one has to solve the more practical problem of how to determine 
action variables in the first place. A very useful tool for that is the Hamilton-Jacobi formalism: Given 
an Hamiltonian H as a function of p and q, the goal is to find a function with 

det a ^ ^ 0 and ^ = Pn such that 

oqj oqj J 

H{q^,-,q^,—= K{h^,-,h^) . (46) 

oqi 

Note, S generates a canonical transformation since 

d5 = d(/j + d/ij) 

3 

^ ^ dpj A dqj = dhj A d0j (47) 

3 

with 0j := Suppose one can find such a solution S to the first order PDE (1461) then, by (1471) . one 

has also found f independent constants of motioifl hj in involution. If det ^ 0 as required then 

dS ^ 

dqj 


hj = /ij(p,q) can be extracted out of the equation ^ = pj due to the inverse function theorem. Apart 
form that also the dynamics of 0,- simplifies, i.e. 


dt ^ dhj 

By setting hf) = E = hi the dynamics becomes especially simple, instead of f constants of motion 

one now has 2f - 1 constants. This seems rather odd, notably if it would hold globally, as it effectively 
reduces the system to a one-dimensional free-particle. Astonishingly, the Hamilton-Jacobi equation (|46l) 


with K = E always has a local solutioij® given by the classical action functional (see e.g. 26-28|). This 
demonstrates that in the above construction the global properties are essential. 0n the other hand, one 
can also fix K(hi,---,hf) = Yhj hj locally which proofs the lemma: 

Lemma 4. For each generic point m e A4 exist a neighborhood Um on whieh the funetions 


Wj := ^Jhj e*^^^ and Wj '■= y/hj ( 


-iujj ©, 


j = 


are well-defined, funetionally independent, canonical conjugated (^{wj,Wfc} = iujjSjk) and obey 

H = ^Wj Wj and {H,Wj} = iujjWj . 

3 

Similarly one can introduce eoordinates Qj '■= ^cosojjQj and Pj := -^J^^sinujjQj which are well- 
defined, funetionally independent and canonical conjugated on Um- In terms of this coordinates Wj can be 
rephrased as 


1 


Wi 


-iC 


\l‘^^3 


Q- with C = - 
2 


® We here use the term ‘constants of motion’ to indicate that they are in general not globally defined. 

® That means that a solution exist in a neighborhood of any generic point m e A4 on which H is not extremal. 
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As already mentioned the Wj’s will in general not define good coordinates on M but can only be 
defined locally. Below it will be demonstrated on the simplest example of a free particle what happens 
if the global properties are ignored (see section IIVB 21) . In contrast to that the complex parametrization 
defined through proper action-angle coordinates is at least well defined in the neighborhood of the whole 


§ 


torus 

In order to understand better what are the necessary criteria for the existence of well-defined 
parametrizations Wj it is a good idea to investigate further the relation between the w-variables and 
those obtained by the complexifier method. The question is whether zj = qj can be obtained 

by lemma |4] given that z define good coordinates, i.e. they are everywhere defined and functionally 
independent, and given that 


{H,zj} = ifj{z) . (48) 

In general, Zj will not be equal to Wj since Wj is canonically conjugated to its complex conjugate Wj but 

{zj,Zk} = {e*^^c {q-^eT^^^xc q^] 

is generically not even constant (see section fill Bl for an example). The momenta conjugated to Zj and 
'Zj are Ilj = e~^^xc(^pj +Uj{qj)) and Yij = e^^xc(^ipj +Uj{qj)) where uj is any function that should only 
depend on qj to ensure that the momenta are Poisson commuting. Note, that the Uj are not arbitrary 
but, due to the stability criterion (|17l) . must be such that 

^ = +5i(z)) = + 5(z)) (49) 

j j 

with 5j(q) = J (?(q) - Uj{qj). In general the functions fj, uj and Qj can depend on complex parameters, 
which is why fj{z) = fj{z). 

Even though the z and w parametrization are generically different they are still closely related. To see 
this let us first derive a set of f constants of motion of (1491) . A function H]^ is a constant of motion of (1491) 
iff 

{e*^xc H,e^^xc H„} = /,(q)p'' + 9(q), = 0 

j 

where := F^xc Hk- Additionally, the maps must be in involution which suggests the ansatz Fi := 
H and Fj = Fj{q) for j = 2,---,f. This yields a system of f - 1 linear, first order partial differential 
equations (PDE) 

^ In fact, this parametrization only breaks down at separatrices where Mh even ceases to be a manifold. Such separatrices 
divide the phase space into several regions on which the level sets may have different properties. An example for that is 
the mathematical pendulum where the phase space divides into an oscillatory and a rotationally branch (see e.g. |^1. 
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that can be solved locally by the method of characteristics (see appendix iBl) . A local solution of such an 
equation depends on a set of f arbitrary initial functions so that it should in general not be problematic to 
find f - 1 functionally independent ones at least for a sufficiently small neighborhood. Yet in practice, it 
can be very hard to actually determine them and it might be easier to use the Hamilton-Jacobi formalism. 
In the following, let 

F(p,q) := H +5j(q)) (51) 

j 

and 


F(p,q) :=e H = f ji^i) p’ + gji<l)) ■ 


(52) 


As before, we want to find maps and S{qi,---,qf,Fi,---,Ff) that satisfy ^ ^ = pj, 


dS dS^ dS dS, 

dqi dq^ j ■’ dqi dq^ j 


(53) 


If det 


^ 0 and det ^ ^ 0 then S and S are generators of canonical transformations that map 


dqjdFfg ' " — Qq^Qp^ 


(p,q) on the conjugated pairs (F,$) and (F, 0) respectively where := and 
variables the reality condition (1491) can be replaced by 


as 

dFi 


> 4 . In these 


e Fj = i fj{z) Hj + gj{z) = -i /^(z) lij + gj{z) = Fj := hj . (54) 

Together with the angles 

©1 - \ ) (55) 

the action variables hj build a new canonical pair (h, 0) because 

{/i„0fc} = ^ [e-^xc{F,-,$J + e*^-c{;p.,4}] = , 

which can be used to construct a w-parametrization as in lemma ID Since and (pj only depend on q 
and F and F the new angle coordinates are of the form 

©j = ^[^i(z,h)+^j(z,h)] . (56) 

It still remains to check whether Wj = \fhj are good coordinates that are everywhere defined and 

independent. For this it is unavoidable to solve the Hamilton-Jacobi equations (1531) . This is particularly 
easy if the Hamilton-Jacobi equation is completely separable, that is, if F is such that F = Y,j FjiQj^ f^)- 
In this case, S = Y,j Sj{qj-,Fj) and S = Y,j Sj{qj-,Fj) with 

Sj{q,F) = Cj{F)+ f 

Jqo 

Sj{q,F) = Cj(Fj)- f 
Jqo 


( 57 ) 
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are solutions of (1531) . Here, Cj and Cj are arbitrary functions which will be set to zero in the following. To 
simplify the notation also the label j will be left away. From ansatz (15 7p one can deduce that 

1 

'qo Jqo 


which yields 


Thus, 


n I - — n 

^{q) =-i / dx——— and <p{q) = ^{q) = i / dx =■ 
Jgo /(*) Jqo /' 


(58) 


0 


-- [ dx^—- [ dx-?— 

2 [Jzo fix) Jzo fix) 


d0: 




(59) 


is well-defined for all points m e on which /(m) = fiq) t 0. Moreover, since the Hamiltonian is real it 
will in most cases be a function of zz so that 


dw A dw = ^d0 A d/i = — ^/(z) ^ -^ + fiz) ^jdzAdz (60) 

is most likely non-vanishing and w = \/he*^® will provide a good parametrization if z does. 

The last thing which remains to check is whether the level sets (14ip are compact, that is, whether (©, h) 
are proper action variables. Note, that if 0 is not a proper angle then w is multivalued and therefore the 
parametrization will break down on the points where 0(m') = 0(m)+n7r. Concluding, if (|59l) and (I60p are 
well-defined nowhere vanishing differential forms then © is most likely a parametrization of an invariant 
torus (for a counter example see section IIV B 21) . 

Even though the above considerations do not prove that integrability and compactness of the level 
sets is a necessary condition for stable coherent states it nevertheless uncovers a close relation between 
existence of (global) action-angle coordinates and those states. In particular, one can only hope to find a 
good global z-parametrization if the system is integrable, has compact level sets and a degenerate dynamics 
(det = 0). 


B. Examples 

To illustrate the construction principle described in section IIV Al and to emphasize its advantages and 
drawbacks, some examples will be discussed. 


1. Radial oscillator 


Apart from the harmonic oscillator the radial oscillator with Hamiltonian 

TT i -2 
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and phase space M = {{p,q) e M?\q t 0} was the only one-dimensional model with stable coherent states 
found in section UlI Bl To obtain stable states one has to consider the complexifier 



for which 


H = iojpq + g{q) ■= F 


and 


2 ;= e q = 






-2 




Starting with that, it is possible to construct a w-parametrization displaying all the properties mentioned 
in lemma |4] by applying the Hamilton-Jacobi method with generating function 


rq I 

= -i dx — {F-g{x)) . 

Jqo WX 


One easily verifies ~ ~u)q ^ ^ all m e AJ, ^ = p and F(q,^) = F. The angles, which are 

canonically conjugated to F and F, are 


, <95 i dSi 

5> = —— =-In X and tP = —= = — In x 

dF LX dF LX 


and the angle conjugated to h is 


0 :=-(Inz - Inz) 

2lx 


This shows that the w-parametrization is given by 


hz 


h 


w:=VhF‘^^ = \r-^ = \ ^z=. ^ 

z y zz ^ v'/i 2 /a ;2 _ 2A 


lx / i /2 


where the last equality follows from 


(61) 


zz = 


\ LX‘^ ) LX^ \LX^ LX^ ) 


\1 

2 

\h^ 1 



^-2A 

/ 

LX 



Of course, the functional dependence of w and w on the original parametrization (p, q) is much more 
complicated than for z and z. However, their algebraic properties are nicer. By construction {w, w} = ilx, 
while for 2 : and 2 holds 




h 


4i 


h 




H 
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Since the Hamiltonian is strictly positive on the phase space also ^ - 2A must be grater than zero. 
Because of that and since z and z are well-defined and independent, w and w are everywhere defined and 
independent. Besides that, the level sets are compact. 

The complex parametrization w can be as well brought into the complexifier form defining 


Q:= 


/i/4 




2A 


Re( 2 ;) and P := 


/i/4 




2A 


Im(2;) 


and C 


= pV 2, Then u, = m Q , _ ip)^ 


Instead of starting with the z-parametrization it is equally well allowed to construct directly the action- 
angle coordinates of the model. With the generating function 

S{q,h)= f dx \/2(h - V(x)) 

Jqo 

the Hamilton-Jacobi approach leads to = h, ^ = P and 


~ dS n dx 1 . y 

9:=-—- / — - = —arcsm— 

dh Jq^ yj2{h-V{x)) 2a; H 


(62) 


where = h? ju:'^ - 2A, y = ojq^ - hjoj and q^ is fixed such that y^ = 1 for simplicity. Since arcsinx = 
-iln(ix + vl-x^), the w-parametrization, which one obtains from this ansatz, is equal to the above up 
to a complex phase depending on the value of q^. Here, w = Yet, the integration of (1621) is more 

involved than that of (|61l) 


2. Free particle 


The free particle with Hamiltonian H = ^ ^ good example for what can go wrong if the global 

properties are ignored. Locally the Hamilton-Jacobi equation, H(q, = h, is solved by the generating 
function 

S{q, h) = yjh q 


with conjugated ‘angle’ 


^ - 1 
dh p 


which is only well-defined for p ^ 0. This yields 

w := = p(cos qjp + i sina;p/g) = 1 Q 

V 2a; 

where Q = ^^pcos-, P = —j^psin-, and C = P^/2. Locally, the differentials dw and dw are well- 
defined and 


dp Adq = dh A d© = dP a dQ = — dw a dw 

ioj 
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is non-degenerate. Bnt it is not possible to extend that to all of M since, firstly, the whole parametrization 
is ill-defined for p = 0 and, secondly, the coordinates (w, w, P, Q) are mnltivalned, that is, w{q,p) = w{q',p) 
for uj- - 27rn + u— and n e N. Even more severe, the motion of w seems to be periodic: 

w(t) = w{po,q{t)) =po[cos{ujt + qolpo) + ism{u;t + qolpo)] . 

All these problems arise from 0 not being a proper angle since the level sets h = const, are not compact. 
On the other hand, for ui « q/p the parametrization is not to bad since np to second order 

P~-p + C>(a;^) and Q =-ojq + . 

Another interesting property of w is that on the level-sets, i.e. for constant momenta, w are the generators 
of the ^-algebra of qnasi-periodic fnnctionj^. The spectrnm of this algebra constitntes the so-called Bohr 
compactification that lies at the heart of Loop Qnantnm Cosmology (see e.g. 29l-l32|L 


3. Anharmonic oscillators 

2 

For a generic one-dimensional Hamiltonian of the form H := ^ + V{q) one can always solve the 
Hamilton-Jacobi eqnation locally by choosing 


S= r dxy/2{h-V{x)) . 

Jqo 


(63) 


To conclnde, we will discnss two examples of an anharmonic oscillator that are widely stndied in the 
literatnre: the Poschl-Teller potential [33 1 


Vp.T.iq) = (cosh(g)) 2 , 

which is an effective potential to describe vibrations in diatomic molecules and for which the Schrodinger 
equation is solvable, and the quartic potential 

Mq) = + y/ , 

which is a standard example in perturbation theory. 

For large q the Poschel-Teller potential is exponentially suppressed since for large q cosh(( 7 ) pa e^^^/2 and 
hence the level sets are compact. Moreover, a local w-parametrization is conceivable. Ansatz (1631) yields 

dS dx coshx 


0 := 


f 

Jo \/2/icosh ^X + A(A + 1) 

,(2hr^n r 

Jo Vsinh^ X + 


This are functions / of the form f{x) = fj for some finite label set T and kj e R. Here, kj = uJjlpo. 
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where 


w ■ 




j^sinh q + \JQ? + (sinhg)^ 


n 


x/27i 


In the last equation the identity sinh”^ q = ln[g + \J\ + g^] was used. The corresponding P, Q-variables are 
even more complicated so that the quantization via the complexifier method can cause severe problems. 

While the motion for the quartic potential is also bounded, the angle 0 is even harder to determine. 
In fact, 


e./’ 

J On 


dx 


iqo \/2h- 

is an incomplete elliptic integral of first kinc]^ that is not solvable in terms of elementary functions. By 
replacing x = -bcoscp and 

2h - u?x^ - = A^(x^ + a^){lP‘ - x^), 

where a? - iP' = (w/A)^ and iPuP + = 2h, this integral can be brought into the so-called Legendre 

canonical form 


0 


i r 

P Jqo 


dx 


a/(x^ + aP){b‘^ - x^) 


k 

= —F{4>, k) + c . 


(64) 


Here, 


c is a constant and 


^sin (p 

F{^,k) = 

Jo 


dC 


f 

Jo 


dip 


'0 \/(l-?^)(l - Jo \/l - k"^ sin^ Ip 

A literal quantization of the corresponding complexifier seems hopeless. On the other hand the quartic 
anharmonic oscillator can be treated, classically as well as quantum mechanically, by a perturbation 
theory. Especially the classical perturbation theory makes heavy use of action-angle coordinates of the 
unperturbed harmonic oscillator. It therefore seems much more promising to study the w-parametrizations 
in this context, however, this would go beyond the scope of this paper. 


DISCUSSION 


It was found that a necessary criterion for the existence of stable coherent states is that the classical 
evolution of the variable z = q depends only on z itself and not on its complex conjugate. However, 

it is not possible to determine other models than the harmonic and radial oscillators from ansatz (|29p that 
allow for such parametrizations. This issue was circumvented by invoking the Hamilton-Jacobi formalism. 
In doing so, one can, indeed, construct local complex coordinates for almost any system that display the 


For more information see e.g. [3j, 13] 


9 
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desired properties and are related to the complexifier formalism. In general, these can, however, not be 
extended globally what was explained exemplarily for the free particle. Another issue arises from the fact 
that the resulting complexifiers and parametrizations are in general not analytic functions of the ‘old’ 
variables (p, q) (see section IIVBI) that potentially causes severe problems when quantizing. 


Let us speculate a bit more about the last point: One essential ingredient for quantization is the choice 
of a polarization. Very roughly, the polarization says how to split the classical phase space Ai of dimension 
2 f into an f-dimensional submanifold V whose elements are then represented as multiplication operators 
on an appropriate Hilbert space. Of course, for Ai = T*C the natural choice of V is the configuration space 
C, which corresponds to the usual Schrodinger quantization. In the case of the harmonic oscillator we 
already mentioned implicitly another possible polarization, namely, that defined by the complex parameter 
a = q - ip. This leads to the Hilbert space S)‘^{C,du) of holomorphic square integrable function on which 
the annihilator d acts by multiplication. The resulting quantum theory is unitary equivalent to the 
Schrodinger quantization and the transformation between the two representations is given by the Segal- 
Bargann transform (HOI) . This close relation between coherent states and complex polarizations is much 
more generic and often used in geometric quantization. 

As was mentioned in section Hi Cl the complexifier approach follows exactly the above line of thoughts: 
It intends to provide a link between an arbitrarily chosen complex Kahler structure and a quantization on 


the corresponding polarizations. This relation was examined in great detail in |36l-l39| focussing especially 


on unitary equivalence of the resulting quantum theories. An interesting result of 


37 


is that for compact 


Lie groups G with phase space T*G only complexifiers of the form C = -^ give rise to a quantization that 
is unitary equivalent to the Schrodinger one with Hilbert space L‘^{G,dx). Here, A is the Laplacian of 
the group and dx is the Haar measure. But formally A can be identified with the square p^ of the original 
momentum operators while the complexifiers obtained for the various examples in section IIV Bl show a 
much more complicated functional dependence on p. Thus, one should not expect that a quantum model 
obtained by directly quantizing the systems in the new w-parametrization is equivalent to the Schrodinger 
representation. Even though it is, of course, tempting to try the first due to the non-analyticity of G, 
which one often encounters. 

Apart from the complicated structure of the complexifier it is also questionable whether the above 
strategy can be applied if the w-parametrization is ill-defined. All of these are interesting questions which 
deserve a further investigation but go beyond the scope of this paper. 


Remember, the original motivation for investigating the stability of coherent states originated from the 
heuristic implementation of the Hamiltonian constraint in LQG by the formal expression ([T|). Yet, the 
formalism developed in section HVl can blow up already for ‘simple’ one-dimensional models. The technical 
obstacles encountered in the finite dimensional examples of this paper are of course even more challenging 
in QFT ?s such as LQG. On the other hand, the formalism developed in this paper in principle directly 
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applies to symmetry reduced models of LQG such as LQC. In fact there are two observations which support 
the assumption that the above formalism is conceivable in LQC: The first hint is the appearance of the 
Bohr compactification, which plays a crucial role in the quantization of LQC, in the stable state system 
of the free particle. The second hint is more vague, namely, the stable parametrizations z^'z of the radial 
oscillator form an sl(2,M) algebra together with the Hamiltonian H. But this algebra also shows up in 
the modified dynamics of deparametrized, homogeneous and isotropic models. 
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Appendix A: Stable complexifier states for the radial oscillator 
1. Classical and Qnantnm Properties of the Radial Oscillator 


Note that the potential of the radial oscillator, 

2 \ -2 
y( q) - 2 ^ ’ 

with mass m, frequency uj and coupling constant A diverges in the limit q 0. Therefore, the classical 
motion is confined to the positive (or negative) axes. The equation of motion can be solved, for example, 
by integrating the law of conversation of energy via a separation of variables. This yields 


q{t) = ± 


E 


mu)^ 




muj 


sin + 4>) 


where 7 = 


.2 2 


s/mu. 


- 2A and the phase (p is determined by the initial condition. 


From the form of the potential, one expects the quantized system to be discrete. In the limit A 0, 
the solutions should approach the solutions of the harmonic oscillator. Due to the barrier at g = 0, only 
the odd solutions, i.e. those solutions which have a knot at zero, are allowed. For ‘big’ A the solutions on 
the positive and negative axis should decouple. This can be checked explicitly by replacing the peak at 


q = 0 through a box potential of width d= 2\ 
a system with energy E is then given by 


and height Vq. The transmission coefficient T for such 


r= 1 + 


Vismh^{^/2{Vo -E)dy 
4E{Vo-E) 


which vanishes for Vq ^ 00 . Therefore, also the quantum system can be restricted to the positive axis. The 
natural choice for a Hilbert space is the space dg) of square intgrable functions with respect to the 

Lebesgue measure dg where q and p are represented as multiplication and derivative operator respectively. 
The corresponding eigenvalue equation. 


H 'tp{x) = E'il){x) 

~ ^ ^ + 2 e j 'ip{x) = 0 


with X := ^ ^ ^ -2 ■ 

solutions are given by a hypergeometric function of the first kind M, that is. 


e := ^ and := 8 mA + 1 is a modified Whittaker equation (see |34| . p. 505) whose 


_ l+Oi 

yj{x) = Ne 2 x 2 M 


1 , , 2 + ck o 

-{2 + a-2e),^-,x 


where is a normalization constant. To ensure that ip is an element of the Hilbert space, the following 
conditions must be satisfied: 
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• lim3;^o^(3j) < oo requires a > -1. Otherwise can not be square-integrable. 

• In order that the Hamiltonian is well-defined, i.e. ||-HV’II < °o, ip needs to be at least once continuously 
differentiable for x ^ 0 (and twice everywhere else). This requires > 1 A > 0. 


Ip is square-integrable iff |(2 + a - 2e) = -re, re e N. Thus the energy spectrum is discrete and 
equidistant 


a. 


: uh{2n + 1 + —) 


(Al) 


In the limit A ^ 0 the eigenvalues approach the odd energy levels of the harmonic oscillator as 
expected and M reduces to an Hermite polynomial (|3j|, p. 505). 


For these specific values the hypergeometric function M becomes proportional to a generalized Laguerre 
polynomial (x), re e N, which build a complete orthogonal system in e“^ dy). Thus, the 

normalized eigenstates are given by 


ipn{x) 


2re! 


\ r(f+ i + «) 


a + 1 

X 2 V, 




(A2) 


Remark: Since the eigenvalue equation for the radial oscillator is a second order partial differential 
equation there exist of course two independent solutions. Here the second linear independent solution is 


ip>{x) = N e 2 x 2 M 


-{2 - a - 2e),-—, 

4V 2 


(A3) 


These functions have to obey analogous normalization conditions: 


Q! 

a < 1 En = ujh{2n + 1 - —) . 

In order that the Hamiltonian is well-defined on these states also requires > 1. Therefore A has to be 
greater than zero. On the other hand, a is equal to ±\/l + 8rreA. Thus either a > 1 or else a < -1. For 
a > 0 we get the solutions above, for a < 0 we obtain ()A3I) . 


2. Coherent States 


Recall, that the complex parametrization z = e *"^^0 q associated to the complexifier 

C=—(lp^+mXq-A 

moj V 2 / 


is given by 


\ a? - mXq~‘^ 

V rnu} 


(A4) 


(A5) 
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where a = The Poisson relations for these variables are 

4 H 

= {H,z} = lioz and {H, z} = -lojz 


zz 


For the operators of the rescaled maps 


^ u: j mJ^ 2 j r -2 

L 3 := —H , L_ := ;= 2 : and L+ := ;= z 

2 ’ 4 V 2 4 V 2 


one recovers the commutation relations of su(l, 1 ): 

[L 3 ,L±] = ±L± and [L-,L+] = L 3 . 


(A 6 ) 


Therefore, it is expected t 
Barut and Giradello (see 


rat the complexifier coherent states coincide with the coherent states defined by 




The easiest way to find the associated coherent states 

ipziq) = e'ft 6 ^{q) 


(A7) 


is to express the convolution 6 in terms of the eigenfunctions of C that are given in terms of Bessel 
functions of the first kind J/ 3 (x) (see e.g. 3^, pp. 358-364). That is, the eigenfunction to the eigenvalue 
hiuc^, 0 0 is 

4>c(x) = ^/CT Jc^/2(cx) 

where x = \/^p- These functions are not square-integrable on M'*' as expected but define a convolution: 


= r 

Jo 


dc4>c(x)4>c(y) . 


(A 8 ) 


Thus, one finds 


i’zix) 


exp 


dcce'zft v^J„/2(cx) J„/2(cy) 

( x'^ + z^\^ ^ ^ 

I-^- yal2{xz) . 


y^z 


(A9) 


Here, I/ 3 (x) = i ^J^(^a;) is the modified Bessel function of the first kind. These states can be also expressed 
in terms of the ei genf unctions of the Hamiltonian by utilizing the relation of Bessel functions and Laguerre 


polynomials (see |34| . p.734). This yields: 


/ , /-s 2 n+^ 

2 00 1 2/\/21 ^ 

i-.M -<=-* S(-ir / ' 

n=o y/n!r(n + q:/2 + 1) 


(AlO) 


For 2 ; ^ z^l\/2 and x ^ x^/2 one recovers the coherent states defined by Barut and Giradello 40| up 
to a constant factor as expected. However, the states flAlOj) are not normalized with respect to the real 
measure dx. Instead, one finds 


IIX2||l 2(R+) - "T e ^ Ia/2 




(All) 
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Properties: The states (lAlOl) are over-complete with respect to the measure 

dfi{z) = [Io^yJ]"^£:K„/ 2 ^yJdrd 6 » (A12) 

expressed in polar coordinates z = re’'^ r > 0, (/> e [“§; f ]■ Here, K,g(x) is a modified Bessel function of the 
second kind. The unusual integration domain of (j) originates from the fact that the analytic continuation 
q z for the radial oscillator is defined only on the upper complex plane. Moreover, the states are truly 
stable under the dynamics, i.e. , 


Since (jAlOp are the eigenvectors of the annihilation operator z, they minimize the uncertainty of 


a 


i2 + (i2)t 
2 


2 ~ 9 

— C + x^ 
00 


(A13) 


and of 


b:= 



— {xp+px). 

UJ 


(A14) 


Appendix B: Partial differential equations of first order 


Partial differential equations (PDE) emerge in all kinds of physical problems and are widely studied 
in mathematics. A detailed treatise would go beyond of this work. It is rather intended to give a rough 
overview of properties and solution techniques for quasi-linear, first order PDE mentioned in the main 
text. For the interested reader we recommend the very detailed book by Hilbert and Courant 2^ or the 
more elementary book by Cohen 411. 

A first order PDE is called quasi-linear if it is of the form 


Y,aj{x,u) 


j=i 


du 

dx^ 


6(x, u) 


(HI) 


where aj and b are continues differentiable functions of x := {xi,---,Xn) and u. The equation is linear if 
aj and b do not depend on u. A solution n(x) to (jBip defines an n-dimensional surfaces, called integral 
surface, whose tangential vectors Vi = ^ satisfy (IBII) at every point with coordinates (x,ti). 

To solve m it is sufficient to determine a family of (n - l)-parametric characteristic curves 
x{s,ti,---tn-i),u{s,ti,---tn-i) that obey 

- = aj{x,u) and — = b(x) 

ds ds 

with initial values 4>j{ti,---tn-i) = Xj{0,ti,---tn-i) and u{0,ti,---tn-i) = For a given set of 

initial data on the (re - l)-dimensional initial manifold C there exist a unique solution of (IBII) iff the 
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functional determinant 


A := det 


ai • 


dx^ 

dXr,. 

dti 

dt\ 

dxi 

dxn 


dtn-l 

does not vanish. If A is zero then solutions of the initial value problem can only exist if C is a characteristic 
manifold, that means, that C is generated by a family of (re - 2)-parametric characteristic curves itself 
(see 28| for details). In this case, there exist infinitely many solutions. 

A generic first order PDE can be always written in the form H(x,p,u) = 0 where pj = Suppose H 
is twice continuously differentiable then solutions can be found by a method similar to the above. Namely, 
one replaces the PDE by the system of ordinary differential equations (ODE) 


dH 

ds 


A dH dxj A dH dpj dH du 


^ dxj ds 


dpn ds du ds 


together with characteristic equations 


dxj dH dre A dH dpj dH dH 

ds dpj ’ ds dpj ’ ds dxj du 


Note, that all these methods are local implying that existence of solutions only holds in an appropriate 
neighborhood of a point where the initial functions and the coefficients in fIBip are well-behaved. 

The general strategy behind the above is to replace the PDE by a system of ODEs which are considered 
easier to integrate. Yet, often, it is exactly the other way around that the PDE is easier to solve, for example 
by the method of separation of variables, than the system of ODEs. The idea of turning a system of ODEs 
into a PDE goes back to Hamilton and Jacobi and will be explained in a bit more detail in appendix ICl in 
the context of classical mechanics. 


Appendix C: Hamilton-Jacobi method and canonical transformations 


This short summary of the Hamilton-Jacobi app roach is mainly based on 26|, more information can 
be found in any good mechanics book, e.g. in |27l |. To keep this small discourse as simple as possible, 
the discussion is restricted to time-independent, one-dimensional systems. Most of the formulas can be 
immediately generalized to models with more degrees of freedom. To also include time-dependent systems 
a bit more work would be required. 

Definition 5 (Canonical transformation). A map g : A4 ^ Ai is a canonical transformation iff 


where g* denotes the pull back and H the symplectic form on Ai. 
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By applying Stokes theorem one can easily show that this condition is equivalent to 


pdq- PdQ = 0 


(Cl) 


for any closed curve 7 in Al. It is also well known that a canonical transformation leaves the equations 
of motion form-invariant, which can be proven by examining the transformation of the one-form := 
pdq- H dt on the extended phase space Ai x . To see why, remember that any alternating two-form in 
odd-dimensions has at least one null-direction, i.e. it has at least one eigenvector to the eigenvalue zero. 
Here, this null-direction is given by the integral curve {p{t),q{t)) for which 

/ / \ / w , 1 dH dq dH dq dH 

0 = d^ll{p{t),q{t)) = dp A dg- —dp A dt - —dg Adt = - —dp a dt) + + —)dq'Adt . 

op oq dt op dt Oq 

This obviously implies the canonical equations of motion and consequently Hi captures the full dynamics. 


Theorem 3 ([26|). Suppose g ■ Ai ^ Ai is a canonical transformation that maps (p,q) to {P,Q) then 
there exist functions K(P,Q) and S{p,q) so that 


pdq - Hdt = PdQ - Kdt + dS , 

— - and ^ ^ 

dt dQ dt dP 

Proof. Since condition (jCll) has to hold for all closed curves 7 the one-form pdq - PdQ is exact, which 
means, that there exist a potential S with dS" = pdq - PdQ. Now set K(P(p,q),Q(p,q)) = H{p,q). This 
proofs the first part of the theorem. The second part follows directly from d^5 = 0. □ 

The function S is called a generating function of the canonical transformation g. For one-dimensional 
models there exist two types of such functioni^: 


A. Suppose det ^ 0 then the momentum can be written as a function of Q and q by the inverse 

function theorem. Inserting p = p{Q,q) in S leads to S{p,q) = Si{Q,q). By comparison of d5 and 
d5i one finds 


dS,iQ,q) dSiiQ,q) 

= p and 


dq 


dQ 


-P . 


B. If det ^ 0 then p = p{P,q). The corresponding generating function S 2 {P,q) is obtained via a 

Legendre transformation, that is, S 2 {P,q) = PQ + S{p,q). A comparison of the differentials yields 

dS2{P,q) , dS2{P,q) ^ 

= p and ---= Q . 


dq 


dP 


10 


For more than one degree of freedom there exist a third type depending on a mixture of Q and P-variables. 
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A given function S 1 /S 2 generates a canonical transformation iff q^qq ^ 0 / §^§p ^ 0- This non-degeneracy 
condition is needed to ensure that Q j P can be extracted as functions of p and q. In more dimensions it 
must be replaced by det ^ 0 or det ^ 0 respectively. 

The simplest example of a generating function is S 2 {P,q) = Pq that gives rise to the identity transfor¬ 
mation. Another application is the Hamilton-Jacobi method. The main idea, hereby, is to transform the 
system such that the dynamics is especially simple. This is, of course, always the case if some coordinates 
are cyclic. Thus, one tries to find functions S'! or S 2 so that 

H{ — ,q) = K{P,t) . 
oq 

Then P is obviously constant and Q{t) = Jg (for specific examples see section HV Bl) . 

This is closely related to the action-angle coordinates that can be introduced for models with compact 
level sets = {{p,q)\H{p,q) = h}. Here, a function 82 ( 1 , q) generating the transformation {p,q) 
is constructed that obeys 

dS 2 {I,q) dS 2 {I,q) , , rrrdS , ,,,, 

-X- =P, -XT-= <?^> and — ,q) = h{I). 

oq 01 q 

For models with one degree of freedom Mh being compact is equivalent with Mh being a closed curve in 
AJ that should be parametrized by (f). This leads to the additional requirements 


I = I{h) and 


d(/> = 27r . 


As shown in 


0 , 


Mh 


the function 


S2{I,q)= [ 

J 'y 


pdq 




where 'yi{qo,q) is a curve in Mh(i) joining qo and q, meets all these requirements. 
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